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In their comment— on our study of phonon heat transport in Nd2Cu04,— Sun and Ando estimate 
that the phonon mean free path at low temperature is roughly half the width of the single crystal 
used in our study, from which they argue that phonon scattering cannot be dominated by sample 
boundaries. Here we show that their use of specific heat data on Nd2Cu04, which contains a large 
magnetic contribution at low temperature that is difficult to reliably extract, leads to an underes- 
timate of the mean free path by a factor two compared to an estimate based on the specific heat 
data of the non-magnetic isostructural analog Pr2Cu04. This removes the apparent contradiction 
raised by Sun and Ando. 

PACS numbers: 72.15.Eb, 74.72.-h 


Recently, we reported the effect of sample size and sur- 
face roughness on the phonon thermal conductivity Kp of 
Nd2Cu04 single crystals down to 50 mK.— At 0.5 K, Kp 
was shown to be proportional to \/A, where A is the 
cross-sectional area of the sample. This clearly demon- 
strates that Kp is dominated by boundary scattering at 
and below 0.5 K. 

In their Comment;^ Sun and Ando (SA) calculated the 
phonon mean free path I from our thermal conductiv- 
ity data, using the formula Hp — C<v>l/5 where 
C = PT'^ is the phonon specific heat and < w > is a 
suitable average of the three acoustic sound velocities. 
Using P — 0.42 mJ/mol for Nd2Cu04, they arrived 
at an estimate of I which is about 50% of the average 
sample width W = l^jApK. Because I <W , they con- 
cluded that phonons in our Nd2Cu04 single crystals did 
not reach the boundary scattering regime at low temper- 
ature. 

Because their criticism is based entirely on a quanti- 
tative estimate of Z, it would seem important to know 
the uncertainty in the parameters /3 and < z; > used to 
arrive at this estimate. SA provide no indication of the 
uncertainty on their numbers. More problematic, how- 
ever, is their use of specific heat data on Nd2Cu04 to 
estimate /3. Indeed, because of the large magnetic con- 
tribution to C (T) at low temperature coming from Nd'^^ 
moments, SA had to make some assumptions to extract 


the phonon component from that low-temperature data. 
An alternate and standard approach is to use the non- 
magnetic isostructural analog Pr2Cu04, for which there 
is no such magnetic contribution. The specific heat of 
Pr2Cu04 below 10 K readily yields /3 = 0.19 mJ/mol 

If we use that value for /3, which is half that used by SA 
for their estimate of ^, we arrive at an estimated phonon 
mean free path which is twice as large. This removes the 
discrepancy which is the basis for their Comment. Note 
also that /3 = 0.19 mJ/mol K'* is in excellent agreement 
with the value of /3 = 0.20 mJ/mol K'' calculated by SA 
using the standard expression based on sound velocities 
(Eq.(2) in their Comment). 

We conclude that as best as one can estimate it, the 
phonon mean free path in Nd2Cu04 at T = 0.5 K works 
out to be roughly equal to the average sample width, so 
that one can expect phonons to be scattered by bound- 
aries below 0.5 K. In as-grown crystals with smooth, 
mirror-like faces, it is then reasonable to expect some 
specular reflection as the phonon wavelength gets longer 
and longer with decreasing temperature. The best way to 
verify this is to deliberately roughen those surfaces and 
see whether the thermal conductivity is reduced. The 
purpose of our article^ was to show that it clearly does. 
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